Correlation Effects on the Temperature Relaxation Rates in Dense Plasmas 
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We present a model for the rate of temperature relaxation between electrons and ions in plas- 
mas. The model includes self-consistently the effects of particle screening, electron degeneracy and 
correlations between electrons and ions. We successfully validate the model over a wide range of 
plasma coupling against molecular-dynamics simulations of classical plasma of like-charged electrons 
and ions. We present calculations of the relaxation rates in dense hydrogen and show that, while 
electron-ion correlation effects are indispensable in classical, like-charged plasmas at any density 
and temperature, quantum diffraction effects prevail over e-i correlation effects in dense hydrogen 
plasmas. 
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I. INTRODUCTION 

Temperature relaxation rates between electrons and 
ions is one of many quantities that must be modeled ac- 
curately in order to predict inertial confinement fusion 
(ICF) [ij . The lack of equihbrium in ICF occurs because 
the fusion alpha particles deposit their energy colhsion- 
ally to the electrons and ions at different rates depend- 
ing on the temperature regime The task is chal- 
lenging because ICF plasmas traverse complex physics 
regimes characterized by collective, quantum and corre- 
lation effects. In practice, there have been attempts to 
splice models from the different regimes d g but they 
are ad-hoc and not validated to the required accuracy. 
Recently, direct many-body simulations [1, 0, 0, IE] have 
been undertaken to validate the various models of relax- 
ation rates. 

Since the seminal works of Landau and Spitzer on 
weakly coupled plasmas 0] , a variety of developments has 
been made to calculate the temperature relaxation rates 
in plasmas. The best established, parameter-free models 
have considered either weakly coupled, non-degenerate 
(ideal) plasmas (e.g. [i3j[lH) or include degeneracy ef- 
fects in the limit of weak electron-ion interactions p^ . 
In spite of recent works [l^ , the effect of particle corre- 
lations on the electron-ion energy exchanges in non-ideal 
plasmas is still not definitely understood 

In this paper, we present a model for the e-i tem- 
perature relaxation rates in plasmas that includes self- 
consistently the effects of screening, electron degeneracy 
and correlations between electrons and ions. The model 
resolves the close and distant particle encounters in a 
self-consistent fashion and does not involve ad-hoc cut- 
offs. We validate the model against molecular-dynamics 
(MD) simulations of classical plasmas over a wide range 
of plasma coupfing. We then apply the model to dense 
hydrogen and discuss the relative importance of degen- 
eracy and correlation effects to the relaxation rate. 



The paper is organized as follows. The model is pre- 
sented in sectionini Our derivation intentionally focusses 
on the temporal evolution of the ionic temperature Ti in 
a two-temperature plasma. Taking advantage that ions 
are classical, an equation of evolution for Ti is obtained 
from momentum integration of the exact kinetic equa- 
tion for the ionic phase-space distribution function. The 
equation obtained (Eq. ^ below) expresses the change in 
Ti in terms of the net work done by the electrons on the 
ions. The latter can be obtained from the ionic and elec- 
tronic density fiuctuations in the plasma. We thus pro- 
pose a model for the density fluctuations that includes 
self-consistently the effects of screening, electron degen- 
eracy and correlations between electrons and ions. Using 
this model, the equation for the ionic temperature be- 
comes a simple rate equation. 



dt 
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where the temperature relaxation rate Vie has the 
Landau-Spitzer form i>ie = uqIhA where vq is a "uni- 
versal" rate and In A is the generalized Coulomb loga- 
rithm that carries the many-body effects taking place in 
the plasma. The model reduces to well-known approx- 
imations in the appropriate limits (e.g. Landau-Spitzer 
formula, Fermi-golden rule formula, Lenard-Balescu for- 
mula.) Several technical details of the derivation of the 
model are presented in the appendices. 

In section IIIII we proceed to validate the predictions 
of the model by comparing them with results of MD 
simulations. Because first-principles simulations of real 
non-equilibrium plasmas are not feasible yet, the vali- 
dation upon truly ab-initio calculations is not possible. 
An approximation used by several authors to cope with 
this consists in performing MD simulations of plasmas 
with semi-classical potentials that mimic quantum ef- 
fects in an approximative way and prevent the unphys- 
ical collapse undergone by classical electron-ion systems 
[E Hi 13 ■ However such semiclassical MD calculations are 
no longer fully ab-initio; besides their adequacy to simu- 
late temperature relaxation is not unquestionably estab- 
lished. The only truly ab-initio simulations of plasmas 
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that can be performed are for hke-charged systems made 
of positively charged electrons and ions immersed in an 
inert, neutralizing background. Highly accurate MD sim- 
ulations of the relaxation rate in like-charged plasmas are 
possible and were recently reported in [8j. We therefore 
apply our model to a plasma of like-charged electrons and 
ions. This is legimate since our formalism is independent 
of the nature of the electrons. Electrons can be quan- 
tum and negatively charged as in a real plasma but can 
also be classical and positively charged: it is only when 
one evaluates the model that quantities pertaining to the 
system studied must be used (e.g. classical vs quantum 
response functions.) As we shall see, like-charged systems 
are as challenging as real plasmas for the theory because, 
contrary to the latter, they exhibits correlation effects in 
all regimes of plasma coupling. They therefore provide a 
good test of the validity of a theory that aims to include 
correlation and strong coupling effects. Moreover, the 
points of differences between like-charged and real plas- 
mas make the comparison very instructive to shed lights 
on correlation effects. For classical like-charged plasmas, 
our model reproduces the accurate MD data and joins 
the weakly to strongly coupled regimes. 

In section |IVl we apply the model to real dense hy- 
drogen plasmas. We shall see that, while e-i correlation 
effects are crucial in classical plasmas at any density and 
temperature, quantum diffraction effects prevail over e-i 
correlation effects in real, dense plasmas. 

Finally, we present in appendix [E] a systematic com- 
parison of our model with the coupled-mode theory of 
Dharma-wardana and Perrot , the Fermi-golden rule 
result [l^, [l^ and with the Lenard-Balescu formula ■ 
In the following, the Fourier transform of a space- and 
time-dependent function a(r,t) is defined as a(k, = 



r, , dr f die 



-if k-r- 



a(r, t), where V is the volume. 



II. THE MODEL 



We consider a spatially uniform, un-magnetized, two- 
temperature plasma consisting of one species of ions 
(mass rrii, number density rii, charge Ze, temperature 
Ti) and free electrons (me, rie = Zui, Zee, Te) in a vol- 
ume V . We assume that, at any time the electronic 
(e) and ionic (i) components of the plasma each may be 
characterized by a temperature Tf,{t) and Ti{t), respec- 
tively. Because of the large difference between me and 
rui , it is indeed reasonable to assume that the energy ex- 
change between electrons and ions occurs on a time scale 
that is much larger that the equilibration times Tq within 
each subsystem a = e,i. 

The (classical) ion dynamics can be completely de- 
scribed by the kinetic equation for the single-particle dis- 
tribution function /i(r, p, t) in the phase space consisting 
of the position r and momentum p , 



dt 



Equation ([T]) is readily obtained as the ensemble aver- 
age (. . .) of the evolution equation for the microscopic 
(Klimontovich) distribution function 

Ni 

7V,(r,p;0 = ^5(r-r,(i))(5(p-Pa(t)), 

a=l 

where ra(i),Pa(0 ^-^^ the position and momentum of 
the a-th ion at time t. Since we are looking at time- 
scales larger than the equilibration times Te_,; within each 
subsystem, (. . .) denotes an average not only over mi- 
croscopic replica of the same macroscopic state but also 
over a time scale of order max(Te^i); the notation 5 A de- 
notes the fluctuations of a quantity A around its average, 
i.e. 5A = A - {A) and {5 A) = 0. In Eq.© /, = {N,) 
and ^F,; is the fluctuating part of the total force acting 
on the ions that is induced by the density fluctuations 
^n^^i ^ J dp 5Nf,^i in the electronic and ionic systems. 
The ionic temperature at time t is given by. 



M(t) - f dppV,(r,p,t), 



(2) 



Using Eqs.([T|) and ([2]), we obtain the equation of evolu- 
tion for Ti, 



dJ\ 
dt 



3n, k f 



(5j,(r,t).5F,e(r,i)) 



(3) 



Here (5j,;(r, i) is the fluctuating part of the ionic cur- 
rent density / dpSNi{r,p,t)p/mi. Without magnetic 
fields, 6ji is longitudinal and is related to Sm through 
the continuity equation: dSni{r,t)/dt = —V • 5ji(r,i); 
In Fourier representation, Sji = ojdni{'k,u!)'k/k'^. The 
term SFie{r,t) is the force induced by the electronic 
density fluctuations 5ne(r,t); in Fourier representation, 
(5Fie(k, (jj) = iliVie (fc)fae (k, uj) where the e-i interaction 
potential is Vie According to Eq.(I2]), the evolution 

of the ionic temperature Ti is determined by the statisti- 
cally averaged work done on the ions by the fluctuating 
force 6Fie exerted by the electrons 

Introducing the expressions for (5ji (k, w) and SFie (k, uj) 
in Eq.Q, we obtain 



dJ\ 
dt 



3nikB V ^ 

k 



(4) 

(ia;a;Uie(^)Ini((5ni(k, Li;)^ne(— k, —uJ)) 



(1) 



where Im denotes the imaginary part. We re- 
mark here that when the system is at equilib- 
rium, Tf. = Ti, the right-hand side of Eq.Q van- 
ishes as it should. Indeed, it follows from time- 
invariance that at equilibrium ((5ni(k, Ci;)(5ne(— k, — tj)) 
is real and equal to (27r)^S'ie(fc, i^) where Sie{k,uj) = 
l/(iVeiV*)^/^/^e*"*(^n^k,i)Jne(-k,0)) is the e-i dy- 
namic structure factor [l5|. 

According to Eq.(l4]), an equation for dTi/dt can be ob- 
tained by modeling the density fluctuations in the two- 
temperature plasma, which we proceed to do as follows. 
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Following Ichimaru [3,|T3|, the fluctuations SNt can gen- need the correlation function for the spontaneous density 
erally be split into two parts, SN^ = SN^'^ + SnI""^\ fluctuations, ((5nL'^(k, w)(5n^'^(k', w')>- Using Eq.® and 
JTV-^*'' represents the spontaneous fluctuations due to the 
discrete nature of the particles and that are present even 
in the absence of interactions. 



Ni 



a=l 



N^^ (r, p; = ^ <5 ( r - - -P- (t - io) ) ^ (P - Pa), (5) 



where and Pa are the particle position and momentum 
at some initial time to- SN^"^"^^ describes the fluctuations 
that are induced by the interactions between the parti- 
cles. Upon momentum integration of Eq.([5]), a similar 
splitting holds for the electronic and ionic density fluctu- 
ations, namely 

Sn^{r,t)^5n^^Hr,t)+Sn^^"''\r,t) {a^e,i) (6) 

where dna"^ (r, t) is the spontaneous fluctuations in the 

particle density of species a and Jria "''^ (r, t) is the den- 
sity fluctuations that develop due to e-e, e-i and i-i in- 
teractions. We shall assume that the fluctuations are 
small quantities, i.e. \Sna\ << 71^. This is a very reason- 
able assumption since each subsystem a is supposed to 
be in "internal" equilibrium at temperature Tq,, and the 
plasma is supposed to be uniform. Accordingly, we ex- 
tend the results of equilibrium linear response theory to 
a two-temperature system and express the density fluc- 
tuations ([5]) as [isl. [l7t 

6no.{k,uj)=6ni'\k,uj) (7) 
+ x'~°\k,uj) J2 Vcfiik) {1 ~ Gap{k,uj)) Sni3{k,Lu). 

In Eq.Q, vll^{k,uj) = Vc,p{k) [l - Gc,p{k,oj)) 5nfi{k,uj) 
is the local effective potential seen by a particle of species 
a due to the density fluctuations 6np, so that the product 
X° (k, Lo)v'^JJ {k, uj), where Xa is the free-particle response 
function of species a, the density fluctuations resulting 
from the interactions of particles of species a with those 
of species /3. When \ — Gai3{k,uj) = 1, Eq.([7]) corresponds 
to the random-phase approximation (RPA), also known 
as the Vlasov or mean-field approximation in plasma 
physics. In this approximation, the effective field seen 
by an electron or an ion in the plasma is the field that 
would be seen by a classical, external test charge embed- 
ded in the plasma. This approximation fails to account 
for the correlations that exist between the particle under 
scrutiny and all the other particles in the plasma. For in- 
stance, it does not account for the fact that the effective 
field seen by a particle must not include the contribution 
from that very same particle. The so-called local field 
corrections (LFC) Gap in Eq.([7|) account for the neglect 
of those correlation (and exchange) effects inherent to the 
RPA approximation. 

The model Eq. ([7]) for the density fiuctuations is used to 
calculate the right-hand-side of Eq.Q. To this end, we 



that at to the positions of two different "spontaneous" 
particles are uncorrelated, we obtain |15i] 

{Snl:\Kc.)6n^;\k\u;')) (8) 
= Scf3i27rfSik + k')S{Lu + uj')n^VS^°^ (fc, to) 

where = ^ J dLe^^t{Sn^^\-]^^o)Sni'\k,t)) 

is the dynamic structure factor of the ideal gas of 
species a at temperature Tq,. The latter can be re- 
lated to the free-particle response function Xa through 
the fiuctuation-dissipation theorem. 

= -—n{ncu/kBTa)lmxl{k,uj) (9) 

with n{x) = 1/1 — for quantum particles and n{x) — 
1/x for classical particles. Introducing Eqs.d?]), dH) and 
(111) in Eq.dH), we obtain 



dt SkBT^n-iV 



E 



du 



Vie{k) 



D{k,uj) 



[1-GUk)] (10) 



xuj [n{huj/kBTi) - n{huj/kBTe)] ImXe(fc, a;)Imx°(fc, w), 
where 

D = (1 - 'UeeXe)(l " UhX°) - UeiUieX°Xi ■ 

and Ma/3 = Vafi{l — Gap) [3]. In obtaining (fTO|) we 
have assumed the static LFC approximation Gapik, lu) = 
Ga[3{k,0), which is enough here given the additional ap- 
proximation performed below; the general result is given 
in appendix 1X1 

Equation (|10p can be further simplified by (a) noting 
that typically rrieTi << m^Te and a = uieTi /rriiTe << 1 
can be used as a small parameter, (b) keeping the lowest 
order term in a and (c) using the f-sum rule to perform 
the u! integral. The details of these operations are given 
in appendix ICl Equation ^TU\i reduces to the rate equa- 
tion 



^ = [T, - Te) , 

where the temperature relaxation rate is 



(11) 



(12) 



1 



poo 

/ dkk-^ 


Vei{k) 


Jo 


ee(fc,0) 



(1-G,,(fc)) 



9Imx°(fc,w) 



duj 



j=0 



Here 6^(^,0) = l-v^^{k){l-Gee{k))xl{k,0) is the static 
electronic dielectric function, Gap{k) = Gapik, 0) are the 
static EEC's. 

Our model (|12p is conveniently rewritten in the 
Landau-Spitzer form 

i^ie = vo In A, 
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i.e as the product of a "universal" rate 
8neZ'^e^^/2TTmeiTii 



vo = 



and of a (dimensionless) generalized Coulomb logarithm 

"(fcBTe)3/2 



InA = 



\l TrrTie 


Tie 


I dkk'^ 


Vei{k) 


10 


ee(fc,0) 



1 

(47rZe2)2 

2 

(1-G,e(fc)) 



9Imx°(fc,c.) 



(13) 



w=0 



Our approach, which consists in modeling the den- 
sity fluctuations in a two-temperature plasma to obtain 
dTi/dt, does not rely on the concept of isolated binary 
collisions and allows us to treat the plasma as a single 
entity and to include self-consistently the electron-ion 
interactions. For instance the collective behavior typi- 
cal of a plasma and in particular the screening of the 
e-i interaction due to both electrons and ions are present 
through the dielectric function D{k,uj) in Eq. pU]) . The 
short-range correlations, which especially affect the con- 
tribution of close encounters, are self-consistently added 
through the local-field corrections (see e.g. the factor 
1 — Gie in the numerator of Egs. fTU)) and[T51) The effect 
of particle statistics, e.g. electron degeneracy, is included 
through the response functions. All these effects are not 
ad hoc constructs but are self-consistently derived. Their 
contributions are analysed in the next sections. 

A detailed comparison of our model with others models 
such as the Fermi-golden rule formula [T5], the Lenard- 
Balescu formula [ij, [2^ and the coupled-mode theory 
of Dharma-wardana and Perrot [l3[ is presented in ap- 
pendix [E] Here we just make the following remarks. 
When only e-i correlations are neglected, Gie = 0, and 
Gee is approximated by its value in the jellium model 
Glf , our model reduces to the so-called Fermi-golden rule 
(FGR) formula obtained within the framework of linear 
response theory assuming weak e-i interactions [l2| . Our 
model can thus be regarded as an extension of the FGR 
formula where the plasma is treated as a whole and e-i in- 
teractions are treated self-consisitently. When the LFC's 
are completely neglected, Gaf3 = 0, our model reduces 
to the result obtained using the Lenard-Balescu kinetic 
equations. 

Finality, it is worthwhile to remark that the Coulomb 
logarithm ([T3|) differs from the Coulomb logarithm enter- 
ing the Ziman formula for the electronic conductivity and 
its extensions to strongly coupled plasmas [l^, [l^, [ISl hi 
that the latters involve, in addition to the terms in the 
integrand of Eg. p^ . the static ion- ion structure factor 
Sii{k). Thus, contrary to the e-i momentum exchanges 
that govern the electronic conductivity, e-i energy ex- 
changes are globally insensitive to the details of the ion 
fluctuation spectrum [T^j (see discussion in appendix [C]) 
This is reminiscent of Bethe's result on stopping power 
[Til . [2l| that the energy loss of a fast charged particle in 
a plasma is fixed by the total number of scatterers. 
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FIG. 1: (color) Coulomb logarithm vs plasma parameter for 
like-charged systems. The red line is a fit to the MD data [^. 



III. VALIDATION OF THE MODEL 

We validate Eq. using ab-initio MD simulations of 
classical e-i systems with a pure Coulomb potential. We 
avoid Coulomb collapse (recombination) by using posi- 
tively (like-) charged electrons and ions immersed in an 
inert, neutralizing background. As explained in the in- 
troduction, this is done because fully quantum mechan- 
ical simulations are not yet feasible. Again, this is legi- 
mate since our theory can also be applied to classical 
like-charged systems. In the classical limit (see appendix 
IB|). we have 



d_ 



Xe(fc,0) = 
and Eq. p2|) becomes. 



(fcBTe)3/2V 2 
Ue/ksTe, 



InA 



dk 1 - G^e{k) 

~k |ee(A:,0)P ' 



(14) 



where ee(fc,0) = 1 -I- {kl,Jk^){l ~ Geeik)) and koe = 
(47rnee^/fcBT'e)^/^ is the inverse electronic Debye length. 

The Coulomb logarithm InA for temperature relax- 
ation in a like-charged, classical plasma [Z^ = Z = 1) 
is shown in Fig[l] from various calculations as a func- 
tion of the plasma coupling parameter g = rL/Xue where 
vl — e^/ksTe is the distance of closest approach (Lan- 
dau length) and A^e — l/koe is the electron Debye 
length. The red line is a fit to results from accurate, 
large-scale MD simulations For H « 1, the MD 
simulations confirm the theories p^ . [llj , which regular- 
ize the divergent collision integrals at small and large 
momentum transfer k, and yield ln(0.765 Adc/^'l) (blue 
line). However, these theories breakdown a.t g > 0.1, 
since they do not satisfactorily describe correlation ef- 
fects and InA < 0. Our model (fT4t. whose results are 
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FIG. 2: (color) 1 — Gie{k) (on a semi-log scale) for classical, 
like-charged hydrogen (lower, full curves) and for real hydro- 
gen with quantum, negative electrons (upper, dashed curves) 
at Tie = 1.6 10^'' cm"^' and for F = 0.01, 0.1, 1. and 10.. 



shown by the black dots in Fig. ^ , not only recovers the 
weak coupling limit hi{0.765XDe/rL) but also is in very 
good agreement with the MD calculations over the whole 
range of coupling. 

Correlations effects are accounted for by the LFC term 
1 — Gie- When correlations are neglected (Ga/3 — 0), 
equation ()14p diverges logarithmically at large momen- 
tum k for all plasma coupling. 



InA = 



dk k^/{k^ 



The integral diverges due to close encounters (large k) 
because the neglect of correlations assumes that the pair- 
distribution functions gafsir) ~ 1 everywhere. However, 
as illustrated in Fig. ([3]), repulsion at small inter-particle 
distances r forces gap{r) to vanish continuously at r = 0. 
In our model, we account for the "hole" in gabi^) by 
the term 1 — Gie{k) since the two quantities are related 
by the Ornstcin-Zernicke (OZ) relation between the pair- 
distribution functions and the direct correlation functions 
Ca(}{k) — —Vai3{k){l — Gag{k))/kBT 15]. In particular, 
the OZ relation implies [2^ 



1-G,e(fc) 



Ze D{k,Q) 

wVzx^(fc,o) 



fc25,e(fc) (15) 



where Sie{k) = ^Jn^i J dk{gie{r) — l)e~*'^ '" is the i-e 
structure factor. We calculate Gapik) and gapir) self- 
consistently by using the hypernetted chain (HNC) clo- 
sure gajsir) = eyip{-Vap/kBT + gap -I- Cafi) in the 
OZ relations. The HNC closure is known to accu- 
rately describe correlations in classical Coulomb systems 
[l^, [23'|. Results for gie{r) and 1 — Gie{k) are shown in 
Figs.® and ^ for various values of coupling. As a con- 
sequence of short-range correlations and for all coupling 



FIG. 3: (color) Electron-ion pair distribution function (on a 
semi-log scale) for like-charged hydrogen (full lines) at F = 
0.001,0.01,0.1,1 and 10. The dashed line shows gie{r) for 
a purely classical plasma with negatively-charged electrons 
at F = 0.01. Here F = e^/aeksT^ = {g/V^f^ is used to 
characterize the plasma coupling., where — (S/Aim^)^^^ is 
the mean interparticle distance. 



parameters, 

g^e{r = 0) = 
and Eq. p^]) implies (see appendix ID|) 



lim 1 

fe^O 



G,eik) = 0. 



(16) 



As a consequence Eq. p^ converges for all plasma cou- 
pling. 



InA = 



dkl~G,e{k) 



Thus, while Debye screening cuts off the integral (|T4)) 
at small k (distant encounters), short-range e-i correla- 
tion effects embodied in 1 — Gie{k) are crucial to provide 
the large momentum cutoff. The results obtained for 
InA (black dots in Fig.(IT])) are in very good agreement 
with the MD calculations over the whole range of cou- 
pling. At small g, Eq.(fT4|) recovers the weak-coupling 
limit ln(0.765A£)e/''L)- At higher couphng, it joins the 
weakly to strongly coupled regimes. We also find (not 
shown here) that InA is insensitive to the ion charge Z 
at constant g, consistent with MD simulations [11 . 

It is remarkable that Eq. (fT4|) recovers the result of Ki- 
hara and Aono (KH) [1^ and rederived recently with 
modern regularization techniques by Brown et al. (BPS) 
[m (see also Q.) Both KH and BPS use sophisticated 
regularization techniques to eliminate the divergences en- 
countered when dealing with Coulomb collisions in a 
plasma with the traditional particle and wave pictures. 
In the particle picture, particles undergo binary Coulomb 



r/TL 

FIG. 4: (color) Electron-ion pair distribution functions (on a 
semi-log scale) for like-charged hydrogen as a function of r/r^ 
where rL is the Landau length for and for coupling 0.001 < 
r = ja^ksT^ < 10.. It is remarkable that, for all F < 0.1, 
gie{r) = 0.5 when r ^ IAvl- 



collisions and the dynamics of each charged species is gov- 
erned by a Boltzmann equation with Rutherford cross 
sections. The collision integrals, expressed as integrals 
over the impact parameter, diverge logarithmically at 
large impact parameter because collective interactions 
between charged particles are not included. In the wave 
picture, those collective (screening) effects are included 
through the plasma dielectric function. In the traditional 
approach of Landau, large momentum transfer are ne- 
glected and particles momenta diffuse in the fluctuations 
of the total electric field of the plasma. As a consequence, 
close collisions are not treated correctly and the collision 
integrals diverge at large momentum transfer. At first 
sight, the agreement between our model for like-charged 
plasmas and [lol [ll| might seem curious since these au- 
thors consider real plasmas with negatively charged elec- 
trons in the so-called classical regime where > Xth, 
where Xth = Tt-/ Vi^eksTe is the thermal de Broglie wave- 
length. This is because the fundamental quantities un- 
derlying KH and BPS theories, namely the Rutherford 
cross section and the dielectric function [l^, [ll[, scale 
like the square of the charges and their results are there- 
fore insensitive to the sign of the electron charge. 

Following Landau and Spitzer's original works Q, we 
note that the result ln(0.765ADe/?'L) can be interpreted 
in terms of maximum and minimum impact parameter 
cutoffs by writing In A = ln(bmax / bmin) with bmax = Xoe 
and bmin = r-L/0.765 « IAvl, respectively, bmax is the 
widespread used maximum cutoff that arises from elec- 
tronic screening while bmin corresponds to the distance 
below which gie{bmin) < 0.5. Figure (H]) indeed shows 
that at small coupling gieir) is always equal to 0.5 when 
r « IAtl. Our approach therefore allows us to under- 
stand the statistical origin of bmin'- close collisions below 
that distance are statistically rare and do not contribute 
to the Coulomb logarithm. 
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FIG. 5: (color) Electron-ion pair distribution functions (on a 
semi-log scale) for like-charged hydrogen (lower, full curves) 
and for real hydrogen plasmas with quantum, negative elec- 
trons (upper, dashed curves) at — 1.6 10^* cmT'^ and for 
r = 0.01, 0.1,1 and 10. 

Finally, our model diverges when applied to fully 
classical and negatively-charged electron; indeed, as il- 
lustrated in Fig[3l since gieir) ex:p{—Vei{r) /kBTe) = 
exp{-ZeZ/rkBTe) ^ oo at r = 0, limfe^oo 1 - Gie{k) = 
— oo at all coupling when Ze = —1., indicating the infinite 
attraction between classical ions and classical, point-like 
electrons. On the contrary, KH and BPS converge in 
the classical limit because the Rutherford cross section is 
independent of the particle distributions, and in particu- 
lar of gie{r). In our model, negatively-charged electrons 
must be treated quantum-mechanically, as shown in the 
next section. 



IV. APPLICATION TO DENSE HYDRODEN 
PLASMAS 

Having validated our model for classical like-charged 
plasmas, we now consider fully ionized (hydrogen) plas- 
mas with negative {Z^ = —1), quantum mechanical elec- 
trons. We define the usual coupling parameter F = 
e^/uekBTe where Oe = {S/inneY^^ is the mean interpar- 
ticle distance, the degeneracy parameter O = ksTe/Ep 
where Ep = h'^kp/2me and kp = (STrne)"'^/'' are the elec- 
tronic Fermi energy and momentum, and — ae/an 
where qb is the Bohr radius. Using the quantum expres- 
sion for ^ImXg(fc,0) (see appendix IB|) . 

|;i«'x°(M) ^ -^^^Ifikm (17) 

where f{k) = 0^/2^^ ^ ^(k' /kl-^./E^)/e^-l ^ ^^^^^y 
adds a Fermi distribution factor /(fc/2) to the classical 
expression used in the previous section. Our model (jl2[) 
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10000 



FIG. 6: Relaxation rate i/ie for dense hydrogen at rie = 1.3e25 
g.cm~^ obtained using (fTS)) (full line), the FGR with Gee ~ 
Gii = (dashed lined), and the widely used Brysk formula 
[3 (dotted hue). 



becomes, 



In A: 



Equation differs from the classical limit ([H]) in two 
major ways. First, it converges even when correlations 
are neglected, i.e. Gaf3 = 0, because /(fc/2) vanishes 
exponentially at large k and cuts off the integral for k of 
the order of 2kF\J\ + 6, 

dfc 1 



A; |ee(fc,0)P 



/(fc/2)<c5o vr,ve 



(19) 



For instance, in the non-degenerate limit >> 1, 
Eq.dH]) gives 



InA = 



(fc2 + 



-r 

ln(0.742V^/r) 



(20) 



This result corresponds to that obtained using the quan- 
tum Lenard-Balescu kinetic equation |25|]. It was more 
recently derived by BPS [III and can also be found in 
[23 |. Following Landau and Spitzer, Eq. ([20|) can be 
rewritten as InA — '^^{bmaxl^min) in terms of the maxi- 
mum and minimum impact parameters hmax — ^De and 
bmin ~ 0.778 At/i, where \th — ^/VrneksT^ is the ther- 
mal de Broglie wavelength; this results is to be compared 
to bmin = 1-4 obtained in the previous section for clas- 
sical like-charged systems. 

In the degenerate limit 9 << 1, /(fc/2) w 6'(2fcF - 
k) and the range of integration is limited to 2fc^ (only 
those electrons near the Fermi surface take part in energy 
exchanges). 



InA = 



2kf 



dk 



1 



In 1 



44 



4fc| 



4fc2 



1,2 



(21) 



For illustration, we used in the last equation the large 
wavelength appoximation ee(fc) ~ 1 + k^p/k'^ where kxF 
is the finite temperature Thomas-Fermi screening length 
[26l |. By comparing Eq.(l2T]) with the Brysk formula [3, 4, 
we remark that Brysk et al. used only the logarithmic 
part. At high density, the second term makes a negative 
correction to the logarithmic part of typically 20 — 30%. 

As shown in appendix [Bl the presence of the converg- 
ing factor /(fc/2) in Ea. pT|) and its absence in Eq. p^ 
are due to the difference in the energy excitations lu — 
k • p/nie + hk^/2me and to — k. ■ p/me, respectively. 
In particular, in the quantum electron gas and for the 
small energy transfers w ~ of interest here, large mo- 
mentum transfers (fc^ >> mefesTe/?!^) characteristic of 
close encounters can only involve electrons in the tail of 
the Fermi distribution. Since the latter vanishes expo- 
nentially at large momentum, the cumulating effects of 
recoil energy and Fermi statistics, hereafter referred to as 
quantum diffraction effect, naturally cuts off Ea. ([T5)l at 
large fc irrespective of the strength of the e-i correlations. 

Second, as illustrated in Figs.([l]) and (O for dense hy- 
drogen, real plasmas exhibit correlation effects that differ 
from like-charged plasmas. Here gieir) varies from 

9^e{r = 0) > 1, 

at the origin to gteir) = 1 at large distance r; moreover, 
9ie{r) — > 1 everywhere as the temperature increases. The 
quantity l-Gie(fc) still obeys the OZ relation Eq^Q [H], 
and varies from 1 at fc = to (see appendix [Pjl 



lim 1 - Gie{k) cx gte{r 



0)>0, 



(22) 



fc \eeikW 



at large fc, as illustrated in Figl^l Thus, in contrast with 
the classical limit, 1 — Gie{k) does not cutoff the integral 
(fT8)) at large fc. The cutoff is instead provided by the 
Fermi distribution /(fc/2). Moreover, the calculations of 
m M, [12] suggest that G,e(fc) < and Gee(fc) > 0, 
and therefore Eq. p^ suggests that correlations tend to 
increase the Coulomb logarithm from its value obtained 
assuming Gap = 0, in agreement with the conclusions of 

We estimated the LFC's needed in Ea. (fT5|) following 
the approach prescribed in [sij for dense hydrogen plas- 
mas. For all the densities studied (lO^'^ < < 10^^ 
cm^'^), a similar behavior illustrated in Fig.(IS]) is ob- 
served. At small temperature, Vie stays nearly constant 
up to fcsTe ~ 0.5Ep at a value slightly higher than when 
correlations are neglected j25|]; indeed, InA scales like 
Te at small Tg , which cancels with the Te mi/Q. Be- 
yond fcsTe > 0.5Ep, the rate decreases and at high tem- 
perature it follows the quantum Lenard-Balescu result 
i^Q In A with In A given by Eq. (|20p , indicating that, in con- 
trast with the classical system studied above, correlations 
do not play any important role. Here, quantum diffrac- 
tion effects play a bigger role than e-i correlation effects in 
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determining , and the FGR formula decently estimates 
the relaxation rates. Note that at these high densities, 
electron degeneracy is always important when the plasma 
coupling is large than unity, i.e. < 1 when F > 1. 
We expect that correlations will play a bigger role when- 
ever gi^{r = 0) (and in turn |1 - G^e{oo)\ oc gie{r = 0)) 
significantly increases while the large momentum cutoff 
2kFV^ + imposed by /(fc/2) also increases; this cer- 
tainly occurs at densities and temperatures low enough 
for bound states to emerge and below which our model 
breaks down [3l| . 

V. CONCLUSION 

In summary, we describe a model for the rate of 
temperature relaxation between electrons and ions that 
treats the effects of electron statistics and particle screen- 
ing and correlations in a self-consistent fashion. Such a 
treatment is necessary in order to calculate the Coulomb 
logarithm without ad-hoc cutoffs and with improved ac- 
curacy for the various physics regimes encountered in 
inertial confinement fusion and stellar interiors. The 
key result in our model for the relaxation rate removes 
the uncertainty in the Coulomb logarithm because it re- 
solves the close and distant particle encounters in a self- 
consistent fashion. The distant encounters are limited by 
the plasma dielectric response for Boltzmann or Fermi 
statistics as needed. The close encounters are limited 
by quantum diffraction effects and short-range particle 
correlations. By treating these effects together and self- 
consistently, the Coulomb logarithm that we obtain for 



a low temperature, oppositely-charged plasma is 50% 
smaller than the often used Brysk Q result and 50% 
larger than the Fermi Golden rule. 

In order to validate our model with ab-initio MD sim- 
ulations, we also applied our fomalism to a plasma of 
like-charged electrons and ions. This is motivated by the 
fact that experiments with real plasmas have not been 
able to provide data of sufficient accuracy to resolve any 
of these issues. The available MD simulations can be 
much more accurate but they are fundamentally classi- 
cal. Moreover, when the interparticle potential is modi- 
fied at short distance to include to some degree quantum 
effects, the simulations are no longer ab-initio. We thus 
abandoned the semi-classical approximation and used the 
real Coulomb force in our MD simulations, but we had 
to make the electron and ion charges alike in order to 
avoid Coulomb collapse. Our model is applicable to such 
a system and our results are in excellent agreement with 
MD simulations. This is a valid test of our formalism 
because the extension to a real and quantum plasma is 
straightforward. 

In the future, we plan to perform MD simulations with 
modfied potentials in order to test the semiclassical ap- 
proximation for both like and opposite charged plasmas 
[32| . In addition, we will compare our model results with 
previous models in ICF experiments. 
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APPENDIX A: GENERAL EXPRESSION 



Eq. fTUl) is obtained assuming the static LFC approxi- 
mation Gab{k,uj) = Gabik,0). In general, we obtain the 
more complicated expression, 



dT, _ 2 



E 



, Vie{k) 



(Al) 



X [We.(fc,C^)5°(fc,C^) -ImA,e(fc,C^)5°e(^,w)] 



with Aa[3 



w«Xi) - UeiUieXeX° and 



Ua/3 = 



0.*), D^{1- UeeX°)(l 

Vapi^ — Gafi) (the star 



denotes the complex conjugate.) Eq. (jAl[) reduces to 
Eq. (fTO|) when Gab{k,uj) = Gab{k,0). In that case, the 
LFC's are real and ImA^^ = Ua/slmXa- To first order 
in the small parameter rngTi/ < rriiTi, both Eq. (jlOp and 
Eg. ljAip lead to the expression (|12p for Ui^. 

Note that the right-hand side of Eq. (|A1[) is equivalent 
to Eq.ll]) and therefore vanishes when = Ti. 



APPENDIX B: IMPORTANT PROPERTIES OF 
RESPONSE FUNCTIONS 

For completeness, we list in this appendix a number 
of basic properties satisfied by quantum and classical re- 
sponse functions. Most of these properties can be found 
(sometimes with different notations) for instance in [isj 
and m. 



Classical and quantum free-particle response 
functions 



We consider an homogeneous system consisting of a 
single species of non-interacting particles of mass m and 
characterized by the particle density n, the (inverse) tem- 
perature P = l/ksT and the chemical potential fi. 



If the system is treated quantum-mechanically, the 
density-density response function x^*'-' (/c, w) is given by 



X(°)(fc,c.) 



dp F{p + hk) - F(p) 



(27r)3 - A£;(p, k) + no- 
where F is the Fermi distribution 

F{p) ^ 



1 -hexp/3(e(p) - ^)' 



(Bl) 



(B2) 



and AE{p, k) = e(p -|- ?ik) — e(p) where e(p) = p^/2m is 
the energy of a particle of momentum p. Equation (|B1[) 
implies 



W°Hfc,c.) = ttJ JP[F{p + nk)-F{p)] 



(27r)3 

xS{huj- AE{p,k)) 



(B3) 



" 2^n4/?A;^'\ + g/3(-^.P+M)' ^^^^ 

with = ^ ± afj- By differentiating (|B4| with re- 
spect to we find 



^inixW(fc,c. = 0) = -n/3^!I|^i/(fc/2), (B5) 



with 



(B6) 



In the classical limit {h 0), Ea. (|B3|) becomes 



(27r)3 

xSihuj - A£;(p,k)) 



(B7) 
(B8) 



where Y = yJmPj^Lo jk, Ai?(p,k) = Tik-p/m and F^-i is 
the Maxwellian distribution 



Fci{v) 



gy^'exp(-/3pV2m) 



(B9) 



The frequency derivative at w = is 



^im;^(o)(fc,^^0) = -n/3^^i. (BIO) 

The only change between the quantum and classical ex- 
pressions for -^Ivny}^^ {k,uj = 0) is in the factor f{k/2). 
The difference stems from the difference in the quantum 
and classical recoil energies Ai?(p,k) appearing in the 
delta functions of Eqs. fBS)) and (|B7)) . viz. 



Ai;;(p,k) = e(p + ak)-e(p) 

= ?ik • p/m + ?i^fc^/2m (quantum) 
= Tik-p/m (classical) 
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In the classical case, energy is always conserved only to 
lowest order in the momentum transfer k. In the quan- 
tum case, the leading term at very large k is ?i^fc^/2m, 
independent of the momentum p. As a consequence, at 
Lo = {), the energy-conserving delta functions in Eqs. (|B3[) 
and (|B7|) brings in the /c-dependent factor /(fc/2) in 
the quantum case and the fc-independent factor Fci(p = 
0) = 1 in the classical case. Hence the additional term 
/(fc/2) in the Coulomb logarithm of real plasmas, often 
attributed in the litterature to electron degeneracy, arises 
from quantum diffraction. 



2. The /-sum rule 

As a consequence of causality, a density-density re- 
sponse function xl^i -2) in a one-component system sat- 
isfies 



X[k,z) = - / duj , Bll 



for any complex number z. Because Imx(fc,ti;) is odd 
with respect to lo, the previous relation implies, 

t ^2n+l 



X(fc,z) = --E^ (B12) 

n=0 

where ujn{k) is the frequency momenty of order n, 

/>oo 

<^n{k) ~ / du! uj" lmxi{k,uj) (B13) 



The first moment lui, which carries information on the 
very short-time dynamics of the system, is independent of 
the interparticle interactions (over very short time scales, 
the particles motion is purely kinetic) and is equal to the 
ideal gas value. 



i(fc) = / diVLjlmx^^Hk,^^) = 



Tmk 



(B14) 



APPENDIX C: DERIVATION OF EQ. (fTT|) 

In this appendix, we show how to obtain Eq. (jlip from 
the more general result Eg. pTI)) . To this end, wc note 
that in most pratical applications a = meTi/miTe << 1 
{nie/mi < 1./1815. << 1 when Te — Ti) and, as a conse- 
quence, the w-integral can be performed analytically by 
exploiting the /-sum rule. Similar but not identical cal- 
culations were performed by Hazak et al. [12] to simplify 
the full FGR formula of Vi,, , and by Boercker et al. [iSl] to 
simplify their extended Ziman formula for the electrical 
conductivity. Througout the appendix we use the general 
properties of response functions recalled in appendix IB] 



First, we rewrite the integrand of Eq. (|T0|) with the help 
of the following quantities, 



Xee{k,Uj) = 



X^e\k,L0) 



^ X^e\k,U;) 

l-u^e{k)x^e\k,uj) ee(fc,w) 



(CI) 



II. \ = xf\k, u:)[l-uUk)x^e\kM] 

XiiK'^-:^) — 



D{k,uj) 



xf\k,u;) ^ 
£i{k,uj) 



(C2) 



where D and defined as in the main text. Xee 

resembles to the response function of the interacting one- 
component electron gas (classical or quantum jellium 
model) while xu is similar to the ion-ion density response 
function of the ions in an electron-ion plasma [15]. Since 
Uab{k) = Vab{k){l — G'a(,(fc)) is real, the imaginary parts 
of Xee and xu are IniXe 



Imxi°V|eeP and. 



ImXi 



Inix, 



(0) 



-I- UeiUielmX ee\Xi' 



(C3) 



Therefore, the quantity ImXglniXi /li*]^ that appears in 
the integrand of Eq. (fTO|) writes 



IniXglmx^ IniXe I^Xi 



(C4) 



ImXeelmXii (/c)Uie(fc)[Imx eej \Xi' 



Because a << 1, ions are much slower than electrons 
and the ionic spectrum of fluctuations lmxii{k, u) resides 
on a very low-frequency range. First, at small fc, where 
collective effects predominate and Y = y/mi(3i/2uj/k << 
1, lmxii{k,uj) rapidly vanishes as the frequency lo ex- 
ceeds the ion plasma frequency ujpi = \J AirUeZe'^ /mi^ 
while lTiiXee{k T Lo) peaks at the electronic plasma fre- 
quency LOpe = \f^irn^e^^Jrne, = ^ m.,lZm(,LOpi » LUp,. 
Therefore the quantities in Eq. (jC4p vanish in this limit 
and do not contribute to the u integral in Eq. (fTTjl . Sec- 
ond, ImxiiikjUj) drops to zero very rapidly as the phase 
velocity uj/k of the excitation exceeds the ion thermal 
velocity kBTi/rrii, i.e. Y = \J mi(3i/2uj /k > 1, since in 
this limit ImxuikjUj) « lmxf'\k,Lu) and 

Imxf ^ ik,u) = -n,p,V^Ye-^" . (C5) 

Overall, the range of frequencies that contributes to the 
uj integral in Eq. (flO)) is characterized by F < 1 and there- 
fore aY << 1. In this limit, Imxee can be replaced by 
its low- frequency limit, 



d 

lmxee{k,uj) ImXee(fc,0) 



(C6) 



d 



,{k,0)\^duj 



lmxi";(fc,0) (C7) 



since lu << ujpe and Imxee(^,w) depends on k and aY 
according to (see Ea. (|B4|) ) 

Imxi°Hfc,c^) = 
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In 



1 + e 



= -ne/3e\/^(ar)e-("^)' 
Accordingly, Ea. (|C4[) becomes 



(quantum) 
(classical) 



Imxi°^ (fc, u;)Imx (fc, uj) 
\Dik,u;)\^ 

d 

« —lmxee{k,0) X Ljlmxii{k,uj) 
ouu 



(C8) 



-U^e{k)Uet{k) ^— ImXee(fc,0)j uo'^ \Xtt{k , 

Direct numerical calculations show that second term 
in the right-hand side of Eq. ljCSP contributes negligibly 
to the Coulomb logarithm (systematically around 0.1% 
of total Coulomb logarithm) and we neglect it from now 
on. 

Keeping only the first term of Eq. ljCSP in Eq. (fTO|) . we 
obtain 



dt 



with 
S{k) 



(C9) 



SksT^Tii y 2kBT( 
dk d 

-—---rVie{k)Uei{k)—lmXee{k,i^ = 0)S{k) 



duLj 



^ksTi J yksTe 
For classical electrons, n{x) — 1/x and 

S{k) = kB{T,-TMk) 

where 



lu\Xti{k,uj) 



(CIO) 



Wl(fc) 



dijjujl-mxii(k, oj) 



(Cll) 



For quantum electrons, hLu/ksTe^i << 1 over the low- 
frequency range outlined before, and therefore 



h(jj 
knT, 



knTp 



flUJ 



(C12) 



and S{k) is again given by Eq. (|C10p . 

The moment LOi[k) can be calculated exactly by noting 
that Xii satisfies the /-sum rule satisified by ordinary 
response functions (see appendix IB 2[) . Indeed, at large 
frequency a; >> 1 and according to Eqs. (|B12p and (|B14[) . 



UJ » 1. 



Substituting Eq. ((CT3l) in Eq.(IC2| implies, 
Xii[k,uj)K 5" w >> 1 



(C13) 



(C14) 



If we assume that xa is causal like its equilibrium coun- 
terpart (see the remark in appendix . then according 
to Eq.dHia), 

Xir{k,u))K. u}i{k)/u'^ uj»\ (C15) 

and therefore, comparing Eqs. (|C14[) and (jClSp . 



uji{k) 



(C16) 



Eq. (jC16|) is valid for both quantum and classical elec- 
trons. Direct numerical evaluation of Eq. ljClip confirms 
the sum rule Eq. (|C16p . 

Finally, for both quantum and classical electrons, we 
find 

S{k)^kB{T,-T,f-^ (C17) 



and 



dTi 

— = -..e (T. - Te) 



(C18) 



where fie is given by Eq. (jl2p 

The high accuracy of the small a = meTi/jruTe ex- 
pansion described here is illustrated in Figs [T] and El for 
classical electrons and for quantum electrons with degen- 
eracy = 1 and 6 = 0.1, respectively. Figs. [7| and [8] 
show a comparison between the fc-integrands in the fol- 
lowing two expressions for the Coulomb logarithm with 
Vie{k) = AnZZee"^ /k'^, namely 



InA = - 



2 {keTef^ m. 



J dk^[l-G,e{k)] J 



duj 



(C19) 

ImXe(fc,w)Imx°(fc,cj) 



\D{k,u;)\^ 

obtained by calculating numerically the uj integral, and 

~{kBnf/^ 



InA = - 



(C20) 



X I dk 
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1-G,e(fc) 



aimxe(fc,w) 



j=0 



ee(fc,0) 

obtained as before by exploiting the f-sum rule. As ex- 
pected, both calculations agree very well. 



APPENDIX D: limfe^oo 1 - G^e{k) 

In this appendix, we prove the results ([TO)) and 
used in the main text. To this end, we combine the OZ 
relation ((T5)) valid for both classical and quantum elec- 
tron with: 

• the asymtotic limit of Sie{k): 
5,e(fc) = ^A^I^ / dr[gie{r)-l]e-'^-'' 



drf{r)e 







47r^neni 



k^ k 



k ^ oo 
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suits of appendix [Bl 




FIG. 7: Dimensionless fc-integrand of Eg. dClQP (full line) and 
Ea. (|C20|) (crosses) for a classical, like-charged plasma ob- 
tained assuming Gap = 0. 



and therefore, using the cusp condition at the ori- 
gin 5-e(0) = —2gie{0) / ub I3l| for real plasmas and 
5-g(0) = for liked-charged plasmas. 



^§fj^.g^e(r-0) (quantum) 



llfgi^){r = 0) (classical) 



the asymptotic behavior of Xe'\k,0): using the re- 



xi°^(o,o) 



(quantum) 
1 (classical) 



oo 



• the asymptotic limit of D{k,0): 
D(fc,0)-1 fc-cx) 

Combining these results in Eq. (jl5p . we obtain 
cx gie(r — 0) > (quantum) 



1-G,e(fc) 



fc ~ OO 



(classical) 



APPENDIX E: COMPARISON WITH OTHER 
MODELS 

In this appendix, we provide a comparison between 
our model and three other models, namely the "Fermi- 
Golden Rule" (FGR) model described by Hazak et al. in 
[l^ and also b y fisl . the Lenard-Balescu (LB) model dis- 
cussed e.g. in 'l4!| and the Dharma-wardana and Pcrrot 
(DWP) model developed in (ISj . 

To this end, the four models are conveniently expressed 
into similar-looking expressions for the evolution of the 
energy density £i — 3nikBTi/2 as 



7f7 p jy^ pj 

—1 = -2h —-r —\v,^{k)fioAN{k,uj)linxT^{k,oj)lrRx'l''''{k,uj) Fermi Golden Rule (El) 

at J {^ny J 27r ' 



— 2^ /|^h.e(fe)|W(fc,..) ,^ I'fS'ilrW^'"^ Lenard-Balescu (E2) 

' (27r)'^ 7 27r ' \1 ~ \vei{k)\^Xe {k,i^)Xi[k,u)\-^ 



^_2n [j% f^lvM'.ANik,.) ^^Xnk,uy)Imxr{k,u;) oharma-wardana-Perrot ^T^ (E3) 

^.2h [^\v.^{k)\'u;AN{k,u;) [1 - G,e(fc)] ImXe(fc, ^)Imx.(fc, c.) 



27r)3y27ri '^^ ' ' ' \1 - \ve,{k)mi - Ge^{k))il - GUk))Xeik,Lu)x^ik,Lu)\^ 

Our model, Eq.dTO]) (E4) 



r 



with AN{k,uj) = n{huj/kBT.i) — n{hu}/kBTe). Indeed, using Imxa = ImXa ~ '>^a{k){l — Gaa)Xa^\'^ 

Equation (|E4|) is nothing but our model Eq.([lOl) writ- in Eq.([lO]) leads to Eq. (lEi|) . RecaU that here Gee and 



ten in terms of the quantities Gu are the static LFC's of a two-temperature, two- 

component electron-ion plasma introduced in the main 
text. 



n ^ - X^°\k,u) 

Xeik-io) = jrr. (E5) 

l-Veeik)il-Geeik))x^^\k,L^) ^ 

The FGR model (|Eip was obtained by calculating the 

x {kuj) = jkjijj) (E6) energy transfers in the first Born approximation in the 

1 - vu{k){l - Gu{k))x'f\k,LL!) e-i interaction. As a result, in (HH), Xe J are the response 
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we remark that 
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FIG. 8: Dimensionless fe-integrand of Ea. (|C19P (full line) and 
Ea. HC20|) (crosses) for a hydrogen plasma at degeneracy O = 
1. (upper part) and = 0.1 (lower part) and F = 1. 



functions of the interacting one-component electron and 
ion plasmas, 



l-Vee{k){l-G°7{k,u;))x^°\k,Lu) 
l~vu{k){l~Grik,u;))xf'\k,u;) 



(E7) 
(E8) 



where G°^p and G°^^ are the LFC's for the one- 
component electron and ion gas at temperature Tg and 
Ti, respectively. 

The LB model (|E2p is obtained from momentum inte- 
gration of the LB kinetic equations The response 
functions in Eq. (jE2|) are 



Xi'{k,u;) ^ 



x\Hk,u) 



X^'\k,u;) 



1 - Vee{k)x'f\k,LL!) 

xf\k,io) 
l-vu{k)xf'\k,c,) 



(E9) 



(ElO) 



The physical interpretation of the three models 
goes as follows. The PGR formula describes the 
energy exchanges between two weakly interacting, 
one-component electron and ion plasmas. The two 
subsystems are independent from each other. For 
instance, electrons do not affect the inter-ionic in- 
teractions; in reality, however, electrons screen the 
i-i interactions and the ionic plasmon excitation be- 
comes an ion-acoustic mode. Conversely, both the 
LB model and our model treat the entire plasma, 
i.e. electrons plus ions, as a single system, in which 
the collisions are due to the interactions via an ef- 
fective, screening potential. In both models the 
effective potential is not an ad hoc construct but is 
derived from a model like (jE4p for the density fluc- 



tuations Sue^i in the plasma. In the LB model, the 
plasma is assumed to be weakly coupled and the 
density fluctuations SueA are accordingly described 
at the level of the random-phase approximation 
(RPA) jl^, which amounts to neglect the LFC's in 
Eq. ([7]) and therefore in Eq. (|E4|) . It is therefore not 
surprising that our model reduces to the LB result 
jHI by setting the LFC's Gap to zero in Eq. ((E4)l . 
In our model, the effect of particle correlations ne- 
glected in the RPA approximation are modeled us- 
ing LFC's in Eq.©. Correlations modify the LB 
model in into two major ways. First, they mod- 
ify the dispersion relation of the collective modes 
in the plasma, i.e. the poles in Eas. (jE2[) and (|E4[) . 
Howvever, as discussed in the main text and below, 
this effect barely affects the relaxation rate when 
meTi/miTe << 1 since in this limit the e-i energy 
exchanges are globally insensitive to the details of 
the fluctuation spectrum. Second, they account for 
the short-range e-i correlations [geiif = 0) > 1 
and therefore 1 — Gie{k — oo) in real plasma, 
3ie(0) = and therefore 1 — Gie(k = oo) = in 
like-charged plasmas) and bring in the contribution 
of close encounters (as discussed in the main text, 
in real plasma, both short-range correlations and 
quantum diffraction determine the effect of close 
encounters.) 

From appendix [Cl it is clear that in the limit 
meTi/niiTe « I both the LB and FGR models 



InA = 

dk 







k |l-i;ee(A;)[l-G(fc,0)]xO(fc,0)| 



(Ell) 
r/(fc/2) 



Comparing our model with the LB and FGR models. 



with G = for LB and G = G°^P(fc,0) for FGR. 
Eq. ljElip was previously derived and validated in 
[l^ for the FGR model and in [H for the LB 
model. As discussed in the main text, our model 
reduces to Eq. (jElip with G = Gee when the e-i cor- 
relation effects are neglected, i.e. G^e = in Eqs. 
(fH)) or p^ . The three models agree when all the 
LFC's are neglected. 
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As discussed in the main text, Eq. (|Elip diverges 
when electrons are treated classically (i.e. when 
/(fc/2) = 1 in Eq. (jEll[ )) because of the inadequate 
treatment of close e-i encounters. Short-range cor- 
relation effects described by the 1 — Gie term in 
Ea.pi|) are crucial in purely classical plasmas at 
any density and temperature to provide the large 
momentum cutoff. The LB and PGR models do 
not describe those correlations properly and ac- 
cordingly diverge at large k for classical electrons. 



Eq. (|Ell[ ) converges when electrons are treated 
quantum mechnically, because f{k/2) vanishes ex- 
ponentially at large k and cuts off the integral for k 
of the order of 2fcirVl + 6- With G = 0, all three 
models lead to the Coulomb logarithm, 



In A 



dk- 



fc3 



— k 



(k^ + 1)' 
ln(0.742V^/r) r << 1 (E12) 



Note that BPS fu] also recover the same result in 
the so-called quantum limit < Xth] see also [2^ . 



The DWP theory (jE3[) was proposed as the non- 
perturbative extension of the PGR theory. DWP at- 
tribute the energy relaxation to the dynamically coupled 
electron and ion modes in the plasma resulting from the 
strong e-i interactions. These modes are accounted for by 
the dynamical response functions xt^ii^, ^) discussed be- 
low. At melting temperature, the DWP relaxation rates 
predicted for metals are an order of magnitude smaller 
than the PGR rates [s^. Por hot dense hydrogen, ac- 
cording to the results recently reported in [7] for densities 
Vs = 0.5 and = 1 and temperatures 50 < T < 4000 eV 
(see their Pig. 2), the DWP relaxation rates are systemat- 
ically lower than the PGR rates, both rates are close (but 
not equal !) at temperatures T > 1 keV and the DWP 
model deviates non-negliglibly and increasingly at lower 
temperatures. These findings differ from the conclusions 
of our model described in the main text. 

The description and the notations in [^, are 
rather confusing concerning the definition of the quanti- 
ties xf^(fcjw) to be used in Ea. (|E3p . 

According to the recent paper Jl, the DWP model 
is simply obtained by including the denominator 1 — 
vfeXT^xT^ in the integrand of the PGP formula 

(see Eq.(9) and related text in 0]), i.e. xt^i — xTi 
Eq. (jE3[) . This seems unlikely since in that case the DWP 
model (IE3[) would be nothing more than the LB model 
(lE2l) corrected by the LPC 0^1% (in practice, 0,[ll[3| 
use static LPC). Then, just like the LB and PGR models, 
DWP would also be well approximated by Eq. (jElip and 
would not significantly differ from PGR. 



According to a previous paper [33] , the xt^i in Eq- (|E3p 
are the electron and ion response functions of a two- 
component plasma, namely 



D 



- ^(0) 



l-z;„„(l-G„„)xi''^ /D 
M 



1 - «ee(l - Gec)xr\ " " Gu)x 

-Ve^V,e{l ~ Ge;){l - G,e)X^°^xf 



With these definitions, it is hard to see how to perform 
analytically the w-integral in Eq. (jE3|) as is possible with 
all the other models using the /-sum rule. Indeed, since 
is sensitive to the low frequency (ion-acoustic) mode, 

the quantity Imxf^Imxf^/|l — \vei\'^xt^x'i^\'^ cannot eas- 
ily be written as the product of two quantities, one that 
resides on a very low-frequency range over which the 
other barely varies. Accordingly, the DWP theory is 
sensitive to the details of the ion fluctuation spectrum 
and find results that differ from PGR in the way recalled 
above. 

Comparing our model with DWP, we remark that: 

• as discussed earlier, both our model and the LB 
model describe self-consistently the electron-ion in- 
teractions in a plasma and therefore the coupled 
modes at the root of DWP theory. The LB model 
describes those modes at the level of the RPA ap- 
proximation while our model brings in correlation 
effects. Again, because in practice electrons are 
much faster than ions, a = nieTi/miTf^ << 1, the 
temperature relaxation rates are insensitive to the 
mode spectrum and both the LB model and our 
model reduce to Eq. (jEll[) . We suspect that the 
prescription of DWP, supposedly based on some of 
the Peymann diagrams discussed in [l^, overesti- 
mates (multiple counts) the effect of screening. 

• unlike our model, the DWP expression (|E3p does 
not include a term like 1 — Gie{k) in the numer- 
ator. As a consequence, when applied to a classi- 
cal like-charged system, the DWP model diverges 
logarithmically at large-k (the integrand scale like 
1/fc at large k as the LB and PGR formulas). The 
DWP is therefore unable to explain the classical 
MD data discussed in section IIIIl while the con- 
cept of coupled-modes at the basis of DWP does 
not only pertain to real plasmas but also applies to 
classical like-charged systems. 

• for dense hydrogen, our model predict rates at most 
50% larger than the PGR, in contrast with the find- 
ings of [7] mentionned above. 



